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Quarkonium suppression in heavy ion collisions is a potential signature of the formation 

■ of the quark-gluon plasma. After a very brief review of the J/ip result at CERN, we restrict 
t^H ' our discussion to the effects of the high-energy multiple scattering of the quark pair in the 

, colliding nuclei. 

: 

§1. Quarkonium in heavy ion collisions 

Investigation of the QCD plasma, which is a system of the strongly interacting 
gas of the quarks and the gluons, is one of the fundamental subjects in the physics 

■ of the strong interaction. 1 ) -4 ) In the experimental study with high-energy nucleus- 
nucleus collisions, the J/ip production is considered as a special signal sensitive to 
the plasma which is created in the early stage of the reactions. The color attraction 
between the heavy cc pair produced in the initial hard parton scatterings will be 
screened in the plasma, which will suppress the binding to the physical J/V>- 

Q_i' This predicted suppression was observed in the NA38 experiment. 6 ) However 

it was seen even for the proton-nucleus (pA) interactions as well as the nucleus- 
nucleus (AB). Because the exponential scaling with the effective target thickness 
L**) is found 7 ) and the 3/ip and the if)' are suppressed in the same rate, 8 ) this 
suppression is interpreted as a normal nuclear effect or absorption of the produced 
pre-resonance state by the nucleons in the nucleus with cr a bs = 6 ~ 7 mb. 7 )' 9 ) The 
large cross section and the state independence are explained as the importance of the 
color-octet pre-meson state in the production of the quarkonia. 10 ) This is a typical 
example to show the necessity of the joint studies of the quarkonium in pp, pA and 
AB reactions. 

The Pb-Pb collisions brought us a real surprise; much stronger suppression was 
reported. n ) Threshold behavior is seen in the data, n )' 2 -* which seems very hard to 
explain theoretically (even with the plasma 12 ) ) . All medium effects must be studied 
duly and critically before the final conclusion is reached, although the plasma is the 
most immediate possibility 

The absorption of the J ftp state by the comoving secondaries produced in the re- 
action is examined extensively 13 ) . As a small color singlet object, J/?/>'s interactions 
with the light secondaries are expected to be weak. Quantitatively, however, it is 
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largely unknown with little experimental information. Model calculations are done 
for the processes like tfm^p — > DD,D*D*, etc. But the resulting cross sections at 
the threshold region are very model-dependent. 14 -* We note that the possible modifi- 
cations of the D state in medium also affect the quarkonium dissociation rate. 15 )' 16 -* 

NA50 collaboration recently reported 2 ) the smaller cross section <7 a b s = 4.4 ±0.5 
mb with the improved statistics, and the J/Y> yield in the peripheral region with the 
smaller errors, which apparently makes the threshold clearer. Any suppression model 
proposed so far should be subjected to this data. 

Strictly speaking, the formation amplitude of the J/ip must be calculated in 
the AB collisions, where the cross section of the J/ip or the pre- meson state with a 
nucleon is unlikely to have a definite meaning. In §2 we study a model in which the 
nuclear effect modifies the cc state before forming the J/ip- 

At the collider energies, the situation will change at least in two respects. One 
cannot think of the independent scatterings of the J ftp off the valence nucleons 
because they are localized in the extremely thin regions in the CM frame (see §3). 
Multiple cc-pair production may alter the assumption of no accidental coalescence, 
which seems valid at the CERN-SPS energy owing to the rare cc production. 17 ) 

§2. A multiple scattering model for 3 ftp suppression 



Recently Qiu, Vary and Zhang (QVZ) pro- 
posed a new nuclear mechanism for the quarko- 
nium suppression, which successfully explains 
the observed data at CERN 1S ) . In their model 
the inclusive cross section of the J/Y> produc- 
tion in the collision of the hadrons A and B is 
written in a factorized form: 

ab— *cc (Q 2 )\ 



a,b \ ^ 

/ dx F 4> a/A (x a )(t) b / B {x b )^2 b ■ 

J X a * 




+ x b 



J/v(<? 2 )> 

(2-1) 
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Fig. 1. 3/tp cross section to 

in the pA and AB collisions n ^ 
compared with the models with 
(2-3) [dashed] and with (2-4) 
[solid], respectively. The filled 
square is the Pb-Pb point. 



where J2 a ,b runs over an parton flavors, Q 2 = 
q 2 + Am 2 , 4> a /A( x a) is the distribution function 
of parton a in hadron A, and x F = x a — x b and x a x b = Q 2 /s. a is the parton cross 
section and Kj/^, is a phenomenological constant. The factor F c5 ^j/ 1 p(q 2 ), which 
describes the transition probability for the cc state of the relative momentum q 2 to 
evolve into a physical J ftp, is parametrized as 



Fc-^j/^q 2 ) = Nj/^ 9{q 2 ) 9(4m" - Am 2 c - q 2 ) [l 



Am 



/2 



4m 2 



(2-2) 



This form includes the effect of the open charm threshold at Am' 2 and simulates the 
gluon radiation effect with the parameter ap > which puts the larger weight to 
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the smaller q 2 . 

The nuclear effect is taken into account as the coherent multiple scattering of 
the pair, which at the leading of the nuclear enhancement may result in shifting of 
the relative momentum in the transition probability, 18 )' 19 ) 

F c5 ^/4f) = F cc^/4l 2 + £2L )> ( 2 - 3 ) 

where L is the effective length of the nuclear medium in the AB collisions. The 
model assumption here is the separation of the multiple scattering regime from the 
later formation stage of the physical resonance. We note here that for a large enough 
L such that q 2 > 4m' 2 — Am 2 the transition probability essentially vanishes due to 
the existence of the open charm threshold (2-2). This apparently gives a stronger 
suppression than the exponential one which follows in the Glauber model. 

The vanishing of the probability for the large L may be altered due to the 
subleading nuclear effects, by which we expect a diffusion of the momentum dis- 
tribution. 20 ) We model the momentum diffusion here by modifying the transition 
probability F c5 _^j/^. The cc pair with relative momentum q produced in a hard par- 
ton collision, will change its momentum to q' after the random multiple scattering, 
and then transforms into the J/tp with the probability F cS ^ i / 1 f ) (q' 2 ). After many 
scatterings, this classical, elementary diffusion process of the momentum results in 
the Gaussian distribution around the initial value q with the variance e 2 L. For 
demonstration, we replace the transition probability by 

F c - c ^/4q 2 ) = {27T J L)3/2 J d'q'e-^F^j/^q' 2 ). (2-4) 

Note that F c5 ^j/^(q 2 ) never vanishes for any q although the average momentum of 

the pair increases as (q' 2 ) = q 2 + 3e 2 L after the multiple scattering. Importantly, 
the transition probability behaves more moderately as F cE ^j /^(q 2 ) °c L~ 3 / 2 for the 
asymptotically large L, which stems from the depletion of the normalization factor. 

In Fig. 1 we show our result on the J/ip suppression calculated using the formula 
(2-1) with the smeared probability (2-4) with e 2 = 0.185 GeV 2 /fm. Other parameters 
are fixed to the same as 18): of = 1 and fj/^ = Kj/^Nj/^ = 0.485. Our model 
reasonably fits the data in the p^4 and AB collisions taken from 11) except the Pb- 
Pb point. The curve bends upward slightly in the semi- log plot. The original QVZ 
model (2-3) with e 2 = 0.25 GeV 2 /fm (dashed line) can explain all the data points 
in Fig. 1. The downward bending of QVZ model is the result of the existence of the 
open charm threshold in (2-2) and the uniform momentum shift (2-3). 

This simple analysis indicates the importance of the subleading effect or the 
momentum diffusion to the L-dependence of the J/ip formation. 

§3. High-energy QQ state passing through random gauge fields 

The coherent multiple scattering of the pair becomes more important at the 
collider energies, where two colliding nuclei are Lorentz-contracted to thin slabs and 
therefore the interaction of the cc pair with the valence nucleons should be treated 
as a simultaneous action. 



4 



H. Fujii 



Let us study a simpler problem of the high-energy heavy QQ pair passing through 
a nucleus, which is modeled as a random filed source, and demonstrate that the pene- 
tration probability of the singlet bound state decays non-exponentially. 21 ) Traveling 
through the target nucleus, the quark-anti-quark pair acquires the eikonal phase, 



U(x,y;A) = W(x;A)W\y;A), 



(3-1) 



accumulated along the path with the transverse positions, x and y, with W(x; A) = 
Vex.p[ig /f^ dx + A a ~ (x, x + )t a ]. This is a matrix in the color space and V indicates 
the path ordering of the product. In the definite 
color basis they are written as 



±Tt(W(x)W\y)), 



U ss (x,y) 



U*(x, y) = 2Tr(t b W{x)t a W\y)). (3-2) 
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Fig. 2. High-energy QQ pair passing 
through a nuclear target. 



Here the color transparency is manifest as 
U ss (x,x) = 1 and U^(x,x) = 0. The survival 

probability of the color singlet QQ state <p{r, z) (r = x — y and z is the longitudinal 
momentum fraction) is written as S ~ | J d 2 rdz(p(r, z)U ss if(r, z)\ 2 . 

As we are interested only in the QQ state, we take a closure of the final target 
states and average the probability over the initial target state. This is performed 
by dividing the path into small zones and averaging over the random gauge configu- 
rations in Eq. (3-1), 21 ) which results in a kernel U(x,y)W(x,y) = K(r,f,R — R) 
with R = |(cc + y) and x the coordinate in the conjugate amplitude. The r (R) 
dependence of K gives rise to the relative (total) momentum diffusion of the QQ 
state caused by this interaction with the random gauge fields. 

In the thin and thick target limits, one finds 
the asymptotic behavior of the survival proba- 
bility, S = J d 2 rd 2 rp (r)p (f)K(r,r,0) with 
Po(r) = J dzip 2 (r,z), as 

(L^O), 



S ~ 1 



L 



S oc 



oo 



]_Lin 

N 2 L 

(3-3) 

The first relation defines Li n . The power- 
law suppression L _1 is the result of the 2- 
dimensional diffusion of the relative momentum, 
which corresponds to L~ 3 / 2 in §2 where we as- 
sumed the 3-dimensional random walk instead. 
The 1/iV 2 is the result of the color diffusion. 

For concreteness we calculated the survival 
probability of a small Gaussian state as shown 
in Fig. 3. Furthermore we speculate an analytic 




Fig. 3. Survival probability of the 
singlet Gaussian state traversing 
SU(3) random fields [cross] is 
compared with Eq. (3-4) [solid] 
and the exponential form [dashed] . 
The first term of (3-4) is shown by 
a dotted line. 
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approximation which is compatible with the asymptotic behavior (3-3): 

2 1 N 2 -2 1 

This formula is found to fit the numerical result quite well, and clearly demonstrates 
the non-exponential behavior of the survival probability. 



§4. Concluding Remarks 



It is a theoretical challenge to predict what will happen at BNL-RHIC before the 
new data show up. 2 )' 4 -* Our elementary demonstration is a example which indicates 
that the J/ip production at BNL-RHIC energy will be different from those observed 
at CERN-SPS energy. At the same time one should think over again the CERN-SPS 
data and the available pA data to construct a unified picture of the quarkonium 
production in the nucleus-nucleus collisions. Such steady efforts will make the study 
with BNL-RHIC and CERN-LHC successful, together with other various observables. 
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